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ABSTRACT 

We present a detailed gravitational lens model of the galaxy cluster RCS2 J232727.6-020437. Due 
to cosmological dimming of cluster members and ICL, its high redshift {z = 0.6986) makes it ideal 
for studying background galaxies. Using ne-w ACS and WFC3/IR HST data, we identify 16 multiple 
images. From MOSFIRE follow up, we identify a strong emission line in the spectrum of one multiple 
image, likely confirming the redshift of that system to z = 2.083. With a highly magnified (/r > 

2) source plane area of ~ 0.7 arcmin^ at z = 7, RCS2 J232727.6-020437 has a leasing efficiency 
comparable to the Hubble Frontier Fields clusters. We discover four highly magnified z ^ 7 candidate 
Lyman-break galaxies behind the cluster, one of which may be multiply-imaged. Correcting for 
magnification, we find that all four candidates are fainter than 0.5L*. One candidate is detected 
at > 10(7 in both Spitzer/IRAC [3.6] and [4.5] channels. A spectroscopic follow-up with MOSFIRE 
does not result in the detection of the Lyman-alpha emission line from any of the four candidates. 

From the MOSFIRE spectra, we place median upper limits on the Lyman-alpha flux of 5 — 14 x 
10“®® erg s“®cm“^ (Sct). 

Subject headings: galaxies: clusters: individual (RCS2 J232727.6-020437) - galaxies: high-redshift - 
gravitational leasing: strong 


1. INTRODUCTION 

It has been decades since the first evidence that clus¬ 
ters of gal axies could be used as cosmic telescopes. |Sou-| 


call (1990) originally proposed the idea to study magni- 
hed background galaxies. Several large programs have 
recently been developed to capitalize on the magnifica¬ 
tion gain. Some of these programs include the Cluster 
Leasing And Supernova surv ey with Hubble (CLA SH, 
PI Postman: HST-GO-I2065, Postman et al. 2012), the 
Grism Lens-Amplified Survey from Space (C LASPS, PI 


Treu: HST-GO-I3459, Schmidt et al. 2014), and the 
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Hubble Frontier Fieldj®^ (HFF, PI Lotz: Program ID 
13495). 

These programs were established largely to better un¬ 
derstand one of the remaining frontiers in cosmology, 
the epoch of reionization, referring to the transition of 
neutral to ionized hydrogen in the intergalactic medium 
(IGM) at redshifts above 6. It is likely that st ar-forming 
galaxies at these redshifts reioni zed the IGM (Robertson 
et al.|20I0[|Bouwens et al.|2 0I2|), though the durat ion ot 
the process is still debatecf (jRobertson et al.||2015|). The 
continued search with the Wide F ield Camera 3 (WFC3) 
on board the Hubble Space Telescope (HST) and with 
future instruments on upcoming telescopes such as the 
James Webb Space Telescope and 30-m class telescopes is 
likely to reveal the exact role of the first galaxies during 
this epoch. 

Here we present a glimpse behind one of the best 
cosmic telescopes, the galaxy cluster RGS2 J232727.6- 
020437 (hereafter RGS2327). RGS2327 is the richest 

and most massive cluster discovered by t he RCS-2 survey 
(Gilbank et al.|20Il jSharon et al. 2015). Original imag- 
ing data from thissurvey revealec . a giant arc (Fig. Q. 


Einstei n radii (iGralla et ^ 201 1|), Sunyaev-Zeb dovTch 
masses (Gralla et al. 201^ Menanteau et al.|20I3|) and a 


strong lensing mass ( Sharon et al.||2015) have since been 


estimated, all corroborati ng the extreme mass of the clus¬ 
ter (M 200 ~ 3 X 10®®M.). Sharon et al. (2015) measured 
the redshift of the cluster to be z = 0.6986 with opti¬ 
cal spectroscopy of several hundred galaxies. The cluster 
is at higher redshift than all of the HFF clusters and 
24 of the 25 CLASH clusters. High redshift cosmic tele¬ 
scopes have a distinct advantage over their lower redshift 
counterparts of the same mass distribution: cosmological 


http: //www.stsci.edu/hst/campaigns/frontier-fields 
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Table 1 

Photometric data summary 


Filter 

texp (s) 

magiin,(lcr) 

HST cycle 

F435W 

4208 

28.7 

15 

F814W 

5626 

28.1 

20 

F098M 

6835 

28.0 

20 

F125W 

6635 

27.4 

20 

F160W 

9247 

27.6 

20 

IRAC [3.6] 

10® 

_ 

_ 

IRAC [4.5] 

10® 

- 

- 


Note. — Photometric Filters used to observe 
RCS2327, their exposure times and Icr limiting 
magnitudes. The limiting magnitudes are ob¬ 
tained from the upper limits of the flux error dis¬ 
tribution found from artificial source simulations. 
IRAC exposure times are approximated from por¬ 
tions of the coverage maps overlapping with the 
WFC3/IR FOV. 


dimming of cluster members and intra-cluster light (ICL) 
decreases the obscured fraction of the telescope FOV. 

Key parameters that influence lensing efficiency in¬ 
clude the mass, geometrical configuration, ellipticity of 
the total mass distribution, amount of substructure in 


the clust er (Meneghetti et a l. 2003) and along the line 
of sight (Wong et al. 2012). An efficient cosmic tele¬ 
scope exhibits large magnification factors (/r > 2) over a 
large area of the cluster. The magnification comes with a 
price: in order to recover intrinsic properties of magnified 
objects, knowledge of the local magnification is required. 

In this paper, we present a gravitational lens model of 
RCS2327 and the first search for z ~ 7 candidates be¬ 
hind it. Sectionj^contains a summary of our new imaging 
and spectroscopic data. In Section!^ we show all mul¬ 
tiple images identified so far in RCS2327, including the 
16 new images presented in this work. We present four 
sub-L* z ~ 7 Lyman-break galaxy (LEG) candidates dis¬ 
covered in our new HST data (Section]^. Finally, we put 
upper limits on the flux and the Lya rest frame equiv¬ 
alent width (EW) of the four dropout candidates from 
our ^ectroscopic follow-up with Keck/MOSFIRE (Sec¬ 
tion . 

Throughout this paper, we assume a flat ACDM cos¬ 
mology with h = 0.7, flm = 0.3, Ha = 0.7, and 
Hq — lOOh km s~^ Mp c~^. All magnitudes refer to the 
AB system ( Oke||1974 l. 


2. OBSERVATIONS AND DATA REDUCTION 
2.1. Imaging data 

The original optical imaging data of RCS2327 were 
taken with HST-I-ACS Wide Field Channel (WFC) in 
F435W (PI: Gladders, HST-GO-10 846). Those data are 
presented in Sharon et al. (2015). We obtained new 
optical and near-IR imaging data with HST-I-ACS and 
WFC3/IR as part of the Spitzer Ultra Faint Survey Pro¬ 
gram (SURFSUP). SURFSUP is a joint Spitzer (PI: 
Bradac, PID: 90009) and HST (PI: Bradac, HST-GO- 
13177) Exploration Science Spitzer program. The obser¬ 
vation and data reduction procedures for our HS T and 
Spitzer data are detailed in Bradac et al. (2014) and Ryan 
etalldMIl). 

Exposure times and limiting magnitudes for the HST 
filters employed in this paper are shown in Table To 


measure HST limiting magnitudes, we mimic the photo¬ 
metric procedure described in Section]^ for our F814W- 
dropout candidates. We first create 5000 simulated point 
sources with m = 25, placing them in random locations 
in a combined E098M -I- F125W -I- E160W detection im¬ 
age. We then run SExtractor (|Bertin & Arnouts||1996 1 
in dual-image mode on each of the HS'F hlters in i'a- 
ble[^ finding that it detects ~4600 of the 5000 simulated 
sources. In the five individual HST filters, where there 
are no sources added, we measure the isophotal flux at 
the locations of each simulated source in the detection 
image. The standard deviation of this flux distribution, 
converted to a magnitude, is the Icr limiting magnitude 
(magJim) quoted in Table We assess the positional 
dependence of the limiting magnitude in each filter. In 
the vicinity of each dropout candidate, the limiting mag¬ 
nitudes differ very little from the global limiting mag¬ 
nitudes reported. We therefore use the global limiting 
magnitudes in Section 

As part of the SURFBUP program, we also obtained 
infrared im aging data from th e Spitzer/IRAG [3.6] and 


[4.5] bands (Eazio et al. 2004|). Due to the larger point- 
spread-function m these bands relative to HST-I-ACS and 
WEC3/IR, a different approach is necessary to conduct 
photometry. As a quick summary o f our IRAC photom¬ 
etry procedure (jBradac et al.||20I4l Huang et al. 2015, 
in preparation) , we use the software package TFIT (Lai- 
dler et al.|2007|) and the segmentation map in WFC3/IR 
FI6UW as the high-resolution prior for mixed-resolution 
photometry. TFIT uses the flux within the isophotal 
aperture of each object as the high-resolution template, 
convolves it with a transformation kernel that matches 
the PSFs in F160W and IRAC bands, and for each ob¬ 
ject solves the flux that optimizes the match to the IRAC 
pixel values. We model all detected objects in FI60W 
within a 20" by 20" box centered at each high-z candi¬ 
date to account for blending in IRAC. In the cases where 
object confusion is too severe to obtain good fits (either 
due to imperfect PSF models or color gradients between 
FI60W and IRAC), we derive the 3cr flux limits for each 
high-z candidate from simulations (Huang et al. 2015, in 
preparation). The simulations are similar to those done 
to measure the HST limiting magnitudes. 

2.2. Ground-Based Spectroscopic Data 
Preexistin g spectroscopic data obtained by [Sharon] 


et al. (2015) confirm the redshifts of multiple image sys¬ 
tems ASzR New spectroscopic observations of RCS2327 
took place on 2013 Dec 15, Dec 17 and Dec 18 (UTC) in 


Infrared Exploration (MOSFIRE, 

McLean et al. 2010 

McLean et al. 

2012) on the Keck 

telescope. (Jbserva- 


three nights, with mostly sub-arcsecond seeing (Tablel^. 

We created a multi-object slit-mask using the MOS- 
FIRE Automatic GUI-based Mask Applicatiorp^ On 
Dec 15 and Dec 17 we used the same mask to observe 
RCS2327. On Dec 18, we used a secondary mask to tar¬ 
get different objects in filler slits. However, we observed 
all dropout candidates on all three nights using the same 
slit configuration. We employed a 2.5" ABBA nod pat¬ 
tern consisting of 180 (120) second individual exposures 


http://www2.keck.hawaii.edu/inst/mosfIre/magma.html 
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Figure 1. Multiple image systems used to constrain the lens model (color-coded by system). The large white circles enclose the four 
primary F814W-dropout candidates. The critical curve (white line) at z = 7 from our best-fit lens model is shown. Also shown are stellar 
mass contours (yellow) for cluster members obtained from Spitzer/IRAC data. The color image is a combination of HST filters: F098M, 
F125W, and F160W. 


in our Y-band (H-band) observations of RCS2327. These 
data were reduced and combined using the publicly avail¬ 
able MOSFIRE data reduction pipeline (DRfQ. The 
reduction pipeline differences, stacks, and rectifies the 
nodded images, creating 2-dimensional signal and inverse 
variance spectra for each individual slit. Spectra from in¬ 
dividual nights were reduced separately due to varying 
observing conditions. 

To extract 1-dimensional spectra from the 2- 
dimensional DRP products, we define a constant (in the 
spatial direction) aperture size. The size is chosen to 
enclose 99.5% of the light smeared by the seeinj^ Pix¬ 
els within the aperture are summed at each point along 
the spectral axis to produce a 1-dimensional spectrum. 

https: //code .google. com/p/most ire/ 

We assume a circular Gaussian profile for the seeing PSF in 
this calculation. 


Flux calibration is done by comparing a 1-dimensional 
extracted spectrum of a standard star we observed to a 
model spectrum of the same spectral type. Model spec¬ 
tra are obtained from the CALSPECFu database. Before 
comparison, we scale the model spectrum to the appar¬ 
ent magnitude of the observed standard, then correct it 
for the airmass and galactic extinction of our observa¬ 
tions. The ratio of the scaled, corrected model spectrum 
to the observed telluric spectrum represents the sensitiv¬ 
ity function of the telescope. 

To obtain calibrated science spectra of each dropout 
candidate, we first extract a 1-dimensional spectrum 
from the 2-dimensional signal spectrum produced by the 
DRP. We then correct this spectrum for airmass and ex¬ 
tinction. The corrected spectrum is multiplied by the 
sensitivity function obtained as described above. We 

http://www.stsci.edu/hst/observatory/crds/calspec.html 
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Table 2 

MOSFIRE spectroscopic data 
summary 


Table 3 

Multiple Images 


Night 

band 

texp (s) 

Seeing (”) 

2013 Dec 15 

Y 

3579 

0.75 


H 

3340 

0.65 

2013 Dec 17 

Y 

3579 

0.88 


H 

3340 

0.92 

2013 Dec 18 

Y 

3579 

1.20 


H 

3340 

0.86 


Note. — Seeing was measured in obser¬ 
vations of alignment stars in J-band directly 
preceding science exposures. 

combine spectra from the three separate nights of ob¬ 
servations to improve the signal-to-noise ratio ( S/N) of 
the spectra . Follo wing an approach described by [FinkeT] 
stein et al. (2013), we scale the noise spectra so that the 
standard deviation of the S/N spectra equals 1. This 
ensures that the noise spectra accurately reflect the er¬ 
rors in the signal spectra. S/N spectra are smoothed to 
the FWHM of an unresolved emis sion line, which is ~ 3 
pixels throughout Y and H-band ( McLean et al.|[^012 ). 

The seeing during our observations is measured in J- 
band alignment images of stars taken prior to each obser¬ 
vation in Y-band. We use the IRAF routine imexamine 
to measure the FWHM of five alignment stars in each 
imag^^ The seeing value for each Y-band observation 
reported in Tablej^is the mean of the FWHM of the light 
distribution measured from five different alignment stars 
in the same frame taken before each set of science expo¬ 
sures. H-band science exposures were always taken di¬ 
rectly after Y-band science exposures. Consequently, we 
calculate the seeing in H-band by comparing the FWHM 
of the continuum trace of a bright object in both bands. 
We account for slit losses due to the seeing after calibrat¬ 
ing our science spectra. 


3. LENS MODEL 

The lens modeling method we emp loy, SWUnited 
(Bradac et al. 2005 Bradac et al. 2009), constrains the 
gravitational potential within a galaxy cluster field via 
minimization. It takes as input a basic initial model 
for the potential. A is then calculated from strong 
and weak gravitational lensing data on an adaptive, pix- 
elated grid over the potential (V'fe) established by the 
initial model. Areas closer to the core(s) of the mass dis¬ 
tribution and in the vicinity of multiple images are re¬ 
constructed with higher resolution. By iteratively solving 
a set of linearized equations satisfying dx^/dtpk = 0 , a 
minimum x^ is found. Derivative lensing quantity maps, 
such as convergence (k), shear ( 7 ) and magnification (/i), 
are produced from the best-fit potential map. 


3.1. Multiple Image Systems 

The strong lensing signal is estimated using positions 
of multiple image systems, which can be identified in im¬ 
ages of the cluster. To determine if objects are multiple 

IRAF is distributed by the National Optical Astronomy Ob¬ 
servatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 


ID 

“J 2000 

<ij 2000 

Z 

F160W 

A.l 

351.87405 

-2.0643629 

2.98® 

23.00 

A.2 

351.87116 

-2.0624877 

2.98® 

23.47 

A.3 

351.85518 

-2.0652714 

CO 

b 

0 

1 + 

0 0 

00 00 

24.07 

B.l 

351.87420 

-2.0723517 

1.42® 

20.60 

B.2 

351.86529 

-2.0669059 

1.42® 

20.63 

B.3 

351.86065 

-2.0682697 

1.42® 

20.77 

C.l 

351.86283 

-2.0645228 

2.083^ 

22.75 

C.2 

351.86062 

-2.0656345 

- 

22.96 

C.3 

351.87950 

-2.0755689 

9 qc+0.31 

24.97 

D.l 

351.86267 

-2.0645370 

9 -1 q + 0.07 
^•J-O_o.o9 

22.75 

D.2'^ 

351.86099 

-2.0653389 

- 

23.63 

D.3 

351.87950 

-2.0754385 

^•-*-^-0.13 

24.16 

E.l'i® 

351.86797 

-2.0731903 

9 qc:+0.15 

^•OD_0.12 

22.61 

E,2def 

351.86785 

-2.0735543 

- 

- 

F.l 

351.87138 

-2.0838888 

-L-uy_o.o7 

22.22 

F.2 

351.86734 

-2.0852578 

-L-Oo_o.o8 

22.53 

F.3 

351.85361 

-2.0772366 

i.56«:i; 

23.37 

0.1==® 

351.86261 

-2.0775600 

- 

24.66 

G.2®® 

351.86165 

-2.0784200 

- 

23.77 

H.l 

351.86199 

-2.0654498 


25.57 

H.2 

351.86025 

-2.0665151 

9 qt:r+0.28 

^•'J'J 0.29 

25.03 


Note. — Column z list the photometric redshift with 68% 
confidence intervals or the spectroscopic redshift, if available. 
Spectroscopic redshift confirmed by [Sharon et al.| | |2015| l. 
Likely spectroscopic redshift from this work. 

Contamination from cluster members/ICL did not permit a 
reliable photo-z measurement. 

May be part of system A. 

® Not used in the lens model. 

^ Detected as same object as E.l in SExtractor segmentation 
map. 


images of the same source, we require that the measured 
photometric redshift (photo-z) of each image falls within 
the 68 % confidence limits of all other images in the sys¬ 
tem. We compute pho to-z distributions for multiple im¬ 
ages using Le Phare ([Arnouts et al.||1999[ jllbert et al. 
2006). We use the Le Phare keyword /_ML, the median 


of the maximum likelihood (ML) distribution, to esti¬ 
mate the redshift of the system. Another requirement 
is that the morphologies of the images are consistent, 
accounting for the parities, distortion and magnification 
predicted by the best-fit lens model. 

Two multiple image systems (A & B) were previously 
identified in imaging data obtained during the RCS-2 
survey and were spectroscopically confirme d to be at red- 


shifts z = 2.98 and z = 1.42, respectively (Sharon et al. 


2015). These two systems provide the strongest con¬ 


straints for our model. We identify 15 multiple images 
belonging to 6 new systems that obey the criteria out¬ 
lined above. We also identify a candidate counter-image 
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Figure 2. 12" x 19" close-up of multiple image systems C and 

D, showing the configuration of the MOSFIRE slit (0.7" x 7") 
used to target image C.l (white rectangle). The color image is a 
combination of HST filters: F098M, F125W, and F160W. 


(A.3) to the giant arc. We measure its pho to-z 


(Table [3]), 
iron et al. 


and w e find that the spec-z measured by Sharon 
(2015) is within the 68% confidence limits. A.3 also pos- 
sesses the expected morphology of a counter-image to the 
giant arc. 

During our MOSFIRE observations of RCS2327, we 
targeted image C.l (Fig. [^. We identify a single H- 
band emission line at 15436 A (Fig. |^. The photo- 
z’s of C.l and C.2 are precise and consistent with each 
other (Fig. 1^. The two best solutions to the single emis¬ 
sion line wavelength based on photo-z are [OIIIJA5007 at 
z = 2.083 or H/3 at z = 2.176, where the uncertainties are 
dominated by the presence of the nearby sky line. The 
line cannot be [OIIIJA4959 at z = 2.112 because we would 
detect [OIII]A5007, the brighter of the [OIII]AA4959,5007 
doublet, at 15581 A with high significance; the line flux 
ratio F([OIII]A5007)/F([OIII]A4959) is always 3 in the 
absence of instrumental or atmospheric effects. The ra¬ 
tio, R = F([OIII]A5007)/F(H/3), for star forming galax¬ 
ies is highly dependent on the gas metallicity. We com¬ 
pute R for a range of metallicities, Q.lZg^ < Z < Zr^, us - 
ing a polynomial fit determined by Nagao et al. (2006). 
We find that for metallicities Z < U.5Zq, K ^ T For 
larger metallicities, 0.5^0 < Z < Zq, the line strengths 
become comparable. If this were the case, then we would 
expect to see both lines in the spectrum. However, such 
large metallicities are improbable at z ~ 2. Therefore, 
it is much more likely that the line is [OIIIJA5007 at 


z = 2.083. 

It is difficult the measure the line flux of the observed 
emission line at 15436 A due to the bright sky line with 
which it is blended. We extract a 1-dimensional S/N 
spectrum of C.l in H-band and show the positions of 
expected emission lines in the case of the bright emis¬ 
sion line being [OIII]A5007. The peak S/N oi the bright 
emission line at 15436 is 5.1. A sky line heavily con¬ 
taminates the spectrum at the location of [OIIIJA4959, 
at which point the peak S/N is 2. However, we do not 
expect a significant detection in the absence of contam¬ 
ination for the reason described above. At the expected 
wavelength of H/3, we see a feature with peak S/N = 3, 
though offset by ~ lA in the rest frame. This feature is 
in a relatively clean part of the spectrum, but does not 


have the convincing negative residuals at the ejected 
locations due to our dithering pattern (Section J2|. We 
expect H/3 to be less or equally significant to [OHTfA4959 
for metallicities in the range O.IZq < Z < Q.5Zq. We 
conclude that the two features at the expected locations 
of H/3 and [OHI]A4959 are consistent with the noise, as 
expected from the predicted line ratios. 

All other emission line scenarios are strongly prohib¬ 
ited by the photometric redshift. However, we examine 
a few other scenarios. If the emission line at 15436 A 
were Ha at z = 1.352, no other strong lines would be 
expected in either band. Similarly, if the line were part 
of the [OHJAA3726-3729 doublet at z = 3.142, no other 
strong lines would be visible in either band. The ob¬ 
served line would have to be the 3729 line, in which case 
the 3726 line would be fully obscured by the sky line. Be¬ 
cause we lack a second line to confirm the spectroscopic 
redshift, we cannot completely rule out either the Ha or 
the [OHJAA3726-3729 scenario. 

Because multiple image systems C and D are lensed 
to approximately the same positions in the image plane, 
they must be at approximately the same redshift. The 
photometric redshift of D.l (Table is consistent with 
the spectroscopic redshift of C.l, as expected. D.2 is 
blended with a cluster member, so we were unable to 
reliably measure its photometric redshift. Systems C and 
D api^ar to be merging, creating a visible color gradient 
(Fig. Using the best-fit lens model, we predict one 
counter-image for each of the two systems. We identify 
two blended counter-images (C.3 and D.3) that display a 
similar color gradient. The photometric redshifts of C.3 
and D.3 are consistent with the spectroscopic redshift 
obtained for C.l. 

Photometric redshift uncertainty is a source of statis¬ 
tic error in our lens model. When available, independent 
measurements of the redshift of multiple images belong¬ 
ing to the same systems can reduce the statistical un¬ 
certainty, and, under certain conditions, mitigate catas¬ 
trophic redshift errors. We do not include the photomet¬ 
ric redshift error in our error budget. Instead, we treat 
the most probable redshift given by the joint photomet¬ 
ric redshift distribution of each system as a spectroscopic 
redshift in our modeling. 

A potentially large source of systematic error in strong 
lensing analyses is the misidentification of multiple im¬ 
ages. While we require multiple conditions to be satis¬ 
fied in order for images to be considered part of the same 
system, we cannot rule out the possibility of misidenti¬ 
fication for systems lacking spectroscopic confirmation. 
Deep spectra needed for redshift confirmation are ex¬ 
pensive and often lacking for the majority of systems 


in strong lensing analyses (e.g. 

Bradac et al.|2008 

Zitrin 

et al. 12009 

Merten et al. 12011 

. 


3.2. Weak Lensing Galaxies 


We generate a weak lensing galaxy shape catal og from 


the ACS F 814W data using t he KSB formalism ([Kaise r 
et al.||1995 ). We largely follow Schrabback et al. 12610 ) 


but additionally apply the pixel-based algorith m tor the 
correc t ion o f charge-transfer inefficiency from |Massey| 


et al. (2014), a revised weak lensing weighting scheme 
based on the measured rms ellipticity as a function of 
magnitude. We remove galaxies with S/N < 10, where 
S/N is defined as the AUTO flux divided by the AUTO 
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Figure 3. Top panel: a portion of the 2-dimensional H-band reduced spectrum of multiple image C.l. A strong emission line is seen at 
Aobs = 15436 A, which is likely [OIIIJA5007 at redshift 2: = 2.0828 di 0.0004. The red horizontal ticks at the left and right edges of the 
spectrum indicate the expected vertical position in the slit from broadband HST data. A skyline can be seen just blue-ward of the emission 
line, limiting the precision of the redshift measurement. The portion of the spectrum shown is 6.7^^ in the vertical direction. Bottom 
panel: 1-dimensional signal-to-noise ratio (S/N) spectrum extracted from the 2-dimensional spectrum. The spectrum is shown at rest 
frame wavelengths assuming redshift 2: = 2.083. Vertical dotted lines indicate the w avelengths of predic ted emission lines in this scenario. 
The line flux ratios are consistent with the line flux ratios for star-forming galaxies (|Nagao et al.j2006|l. Shown below the spectrum is an 
arbitrarily-scaled pure noise spectrum (green). The peaks of the noise spectrum are due to strong sky emission lines contaminating our 
spectrum. The S/N spectrum was smoothed using a gaussian kernel of cr = 1 MOSFIRE pixel. 





Figure 4. Photometric redshift distributions for images C.l and 
C.2. Together, the distributions put a tight constraint on the red¬ 
shift of the system, in excellent agreement with the spectroscopic 
redshift (dotted black line) suggested by the emission line detected 
in the MOSFIRE H-band spectrum of C.l. Note that there are no 
other probability peaks outside of the range shown. 


flux error measured in the SExtractor shape catalog. See 
Schrabback et al. (2015, in preparation) for further de¬ 
tails. 

Red cluster members are removed from the shape cat¬ 
alog by identifying the red sequence. The identification 
is done in the F814W — F098M vs. F098M plane. We 
isolate the red sequence via 

20 < Togs < 27 

1.54 — 0.05 • logs < ^8i4 ~ logs < 1-75 — 0.05 • logsj 


where logs and are ISO magnitudes in the F098M 
and F814W HST filters, respectively. ISO magnitudes 
were measured then scaled to the F160W AUTO magni¬ 
tude (see Section |3. 

After removal oithe red sequence members, 715 galax¬ 
ies r ema in within the leasing field we consider (Sec¬ 
tion 3.31. Each galaxy is weighted by 1/rmSg, where 
rmSe IS the rms of the measured ellipticity. Galaxies are 
assigned a photometric redshift equal to the median pho¬ 
tometric redshi ft, z = 1.18, of the selec ted weak leasing 
galaxy sample. Applegate et al. (20141 show that using 
photo -0 point estimates for weak leasing galaxies leads to 
large biases in the inferred cluster masses. However, the 
signal normalization in our lens modeling code is dom¬ 
inated by strong leasing. For the same reason, we are 
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not concerned with removing the residual blue cluster 
members that normally dilute the weak lensing signal. 

3.3. Magnification Maps 

We construct magnification maps of RCS2327 in 
a 4' X 4' field centered at the location of the 
brightest cluster galaxy (BCG; aj 2 ooo, ^J 2000 = 

23*'27“8".646,-02°04'37.86). The lensing field size is 
chosen to be larger than the FOV of ACS WFC, the 
widest of of our HST observations, to avoid edge effects 
in the mass reconstruction. We provide the lens mod¬ 
eling code with an initial model of non-singular isother¬ 
mal ellipsoids (NIEs) centered at the light peaks of the 
BCGs. The input NIE parameters are velocity disper¬ 
sion, core radius, ellipticity and orientation. We first 
attempt to model the cluster using only two NIEs in our 
initial model. However, we are unable to reproduce all 
of the strong lensing data without the addition of more 
cluster member halos. We find the best fit to the data 
starting from an initial model of five NIEs. We deter¬ 
mine the best-fit run by the strong lensing source plane 
x|l. In the best-fit run, we find XsL/^images = 2.2, 
where Nimages is the total number of multiple images 
input to the model. Magnification maps are obtained 
from the best-fit reconstruction of the potential, which 
is smoothed using a circular Gaussian kernel before ex¬ 
tracting magnification information. 

We compare the magnification properties of RCS2327 
to a sample of the most powerful known lenses, the HFF 
clusters. Eig.f^ shows the cumulative source plane area 
at z = 7 as a function of the magnification for RCS2327 
and the six HFF clusters. The image plane area (~ 4.5 
arcmin^) that is lensed back to the source plane for this 
calculation is approximately the size of the WFC3/IR 
FOV and is the same for each cluster. RCS2327 is most 
similar in magnification properties to MACS0717, yet 
with rounder critical curves which makes it more suitable 
for observations with WFC3/IR. For values of magnifi¬ 
cation /r > 4, RCS2327 magnifies a larger area than 5/6 
of the HFF clusters. The Frontier Fields magnification 
maps we use to make the comparison are created from 
our version I models submitted as part of a special STScI 
call to map the Frontier Fields using ancillary datcp^ 

Magnification error maps are obtained by bootstrap re¬ 
sampling the weak lensing catalog. For each resampling, 
the same initial model and run parameters are used. This 
is done 100 times to improve statistics. Confidence limits 
were obtained by taking the central 68 of the 100 sample 
points. We note that the error maps that are generated 
only account for statistical uncertainties. 


3.4. Stellar mass map 

The Spitzer/IRAC 3.6 p,m image samples close to rest- 
frame AT-band for the cluster, so we use the 3.6 fim fluxes 
from cluster members to approximate the cluster stellar 
mass distribution. We first select the red sequence clus¬ 
ter members brighter than the 25th mag in F814W from 
the color-magnit ude and color-color d iagrams following 
the procedure in|Richard et al.| (|2014|). We select a to¬ 
tal of 311 bright cluster members tor their stellar mass 
distribution. 


http://www.stsci.edu/hst/campaigns/frontier-fields/ 
Lensing-Models 


To create an image with 3.6 pm flux from cluster mem¬ 
bers only, we first create a mask with value 1 for pixels 
that belong to cluster members in the F160W image and 
0 otherwise. We then convolve the mask with the 3.6 pm 
PSF to match the IRAC angular resolution, set the pixels 
below 10% of the peak value to zero, and resample the 
mask onto the IRAC pixel grid. We obtain the 3.6 pm 
map of cluster members by setting all IRAC pixels not 
belonging to cluster members to zero and smooth the 
final stellar mass map with a two-pixel wide Gaussian 
kernel. 


4. DROPOUT SELECTION 


The F814W-dropout search area is limited by the size 
of the WFG3/IR FOV, which is 136" x 123", compared 
to the 240" X 240" grid used in the lens model. To 
conduct the search for F814W -dropouts, we clos e ly fol¬ 


low the met hod d escribed by Goe et al. (2006), Stark 
et al. (2009), and |Hall et ^ ( |2012[ ). Specifically, we 


run SExtractor in dual image mode, using a combined 
F098M -I- F125W -I- F160W infrared image as the detec¬ 
tion image. Eor each object detected in the combined 
infrared image, we measure ISO and AUTO magnitudes 
in each of the five HST filter images (Table [^. Huang 
et al. (2015, in preparation) show that ISO colors per¬ 
form best over the magnitude range in which we are 
interested (25 < m < 28). AUTO magnitudes, how¬ 
ever, are more representative of the total magnitude. 
We therefore use ISO colors, but add to the ISO magni¬ 
tude in each filter the constant term: MAG_ISO_F160W 
- MAG_AUTO_F160W to recover the total magnitude 
in each filter. Magnitude errors reported from SExtrac¬ 
tor are scaled by a constant factor in each filter which 
is derived by Huang et al. (2015, in prepa ration). The 
metho d is similar to what is described by [Trenti et al.| 
(2011) to correct photometric errors in the presence of 
correlated noise. 

Here we describe the procedure for selecting F814W- 
dropouts, i.e., galaxies that show no significant flux in 
F814W or bluer bands. We employ two different color 
and S/N selections and then compare the resulting two 
samples. Our fir st selection meth o d is b ased on the color 


criteri a used by Bouwens et al. (2011) and Hall et al. 
(2012), where they targeted zgso dropouts. I'heir cuts 
are also appropriate for efficiently eliminating contami¬ 
nants from our selection, while simultaneously probing 
the region of color-color space occupied by z ^ 7 star- 
burst galaxies (Fig. |^. The first set of color criteria are: 


r F814W - F125W > 1.0 -f 0.4(F125W - F160W) 

J F814W-F125W > 1.0 
iF125W-F160W<0.5 

( 1 ) 

We visually inspect all candidates that fulfill the above 
color criteria. Objects that are clearly detected in 
F435W and F814W are discarded. We settle on 5 objects 
which show no significant flux above the background in 
both AGS filters. 

Our second selection method is similar, but it uses au¬ 
tomatic S/N criteria in place of visual inspection. The 
color and S/N criteria for this selection are: 
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Figure 5. Cumulative source plane area versus magnification at 2: = 7. RCS2327 is compared to our submitted models of the Frontier 
Fields. The image plane area selected for the comparison is the same for all clusters. For each cluster, we select a ~ 4.5 arcmin^ field 
centered at the center of the F160W pointing. For large values of magnification (tt > 4), RCS2327 magnifies a larger area than 5/6 of the 
HFF clusters. RCS2327 and MACS0717, one of the most massive galaxy clusters known ([Limousin et al.|2012^, follow similar curves. 


' (F814W-F125W) > 1.1 
(F814W - F125W) > (0.8 + 1.1(F125W - F160W)) 
(F814W - F125W) > (0.4 + 1.1(F098M - F125W)) 
(F098M - F125W) < 1.25 
(F125W - F160W) < 0.4 
FSIAW > 25 
S/N < 2.0 in F435W 
S/N < 2.0 in F814W 
S/N < 5.0 in F125W 
, S/N < 5.0 in F160W 

( 2 ) 

7 objects fulfill these criteria, 4 of which were inde¬ 
pendently identified by the first selection method. We 
will refer to the four objects that were selected using 
both methods as the primary dropout candidates. To be 
more inclusive in our spectroscopic sample, we include 
the four additional candidates that were selected by one 
method but not by both. We will refer to these four 
as the secondary dropout candidates. We show cutouts 
of the four primary dropout candidates in Fig. After 
the spectroscopic data were taken, the photometric er¬ 
rors were rescaled as described by Huang et al. (2015, 
in preparation). This does not significantly change the 
photometry of the dropouts. However, we note that in 
Fig.[§ we show the updated photometry for the dropout 


candidates, while showing the original color criteria from 
the first selection method to select the spectroscopic tar¬ 
gets. 

Photometric redshifts were not used as part of the 
dropout selection criteria. They were obtained after the 
selection using Le Phare (Table |^. The full photomet¬ 
ric redshift distributions are shown in Fig. All four 
dropouts favor z ^ 7 solutions, regardless of whether we 
include IRAC [3.6] and [4.5] data. We examine the pho¬ 
tometric redshift distributions more closely in Section 

Fig.[T] shows the location of the critical curve for a 
source at z = 7. 3/4 of our dropout candidates are lo¬ 
cated outside of the critical curve. RCS2327 - 1412 is 
located inside the critical curve at its photometric red¬ 
shift, z = 6.8. We check for an image on the opposite 
side of the critical curve at the expected location pre¬ 
dicted by the model, but we do not detect it in the HST 
data. Because we lack strong lensing constraints near the 
location of RCS2327 - 1412, our model is more uncertain 
in this area. Therefore, the model may not be able to 
accurately predict counter-images at this position. An 
additional counter-image is predicted on the other side 
of the cluster, but would be too faint to detect in our 
current HST data. 

As a result of their large magnifications, all four of 
the primary dropout candidates are intrinsically fainter 
than L* at z = 6.^^ We measure their ages and intrinsic 
stellar masses. Errors in these quantities are obtained by 

We use the ch aracteristic luminosity, L*, determined by 
[Bouwens et al.|2011| 
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Figure 6. Four primary (black squares) and four secondary (black circles) F814W-dropout candidates. The grey shaded region is the 
dropout selection window defined in 0' Because all dropout candidates were not detected in F814W, vertical errors are shown as limits 
(black arrows). The color tracks for z > 6 starburst galaxies including varying amounts of dust are shown in blue, cyan, magenta 
and yellow. The fact that all but one dropout candidate falls near the 2 : = 6.0 points on the starburst color tracks is an indication 
of the limited dynamic range in F814W — F125W color. The outlying dropout candidate is RCS2327-1282. Also shown are possible 
conta minants. Colo r tracks for elliptical (E, black), irregular (Im, red), and Sbc (teal) galaxies wer e obtained fro m the galaxy libraries in 
BPZ | |Bemtez]|2000[ ). Colors for O-M stars (green points) were obtained from the Pickles library (|PicklesJ|1998|. Also shown is a single 
density contour enclosing the portion of the diagram occupied by L,T, and M dwarfs obtained from the A. Burgasser SpeX compilation 
(http: //pono .ucsd. edu/-adain/browndwarf s/spexpr ism/). 


resampling the photometry from a gaussian of location 
and scale equal to the measured magnitude and magni¬ 
tude error, respectively. We do this 1000 times and rerun 
Le Phare on each resample. 

5. MOSFIRE NIR SPECTROSCOPIC FOLLOW-UP OF 
F814W-DROPOUT CANDIDATES 

We observed all 8 F814W-dropout candidates (primary 
and secondary) with MOSFIRE Y and H-band for 3 
hours. We performed a visual search for emission lines 
in all reduced 2-dimensional spectra. We subject each 
potential emission line we visually identify to a series of 
tests. 

1. The line is not significantly offset with respect to 
the continuum position. 

2. The line possesses two negative residuals at the po¬ 
sitions expected from the nod pattern, each with ~ 
half the strength of the line itself. 

3. The line FWHM is larger than or equal to the see¬ 
ing FWHM in the spatial direction and is larger 


than or equal to the instrumental FWHM along 
the spectral axis. 

4. The line flux is comparable in data from individual 
nights. 

5. The line is not due to a known sky emission line. 

6 . The redshift suggested by the wavelength of the 
line falls within the 95% confidence limits of the 
photometric redshift distribution calculated using 
Le Phare. 

After visually inspecting the 2D spectra of all primary 
dropout candidates, we do not detect any viable emis¬ 
sion lines. We do, however, detect a line feature at 3.6cr 
in the combined Y-band data of object RCS2327 - 843, a 
secondary dropout candidate (Fig.|^. The feature is cen¬ 
tered at 10245 A, which, if Lya, would mean the galaxy is 
at z = 7.42. Despite passing tests 1, 3 and 5, the line fails 
tests 2, 4 and 6. Upon inspecting the reduced spectra of 
this object from individual nights, we only see a signifi¬ 
cant detection in the last night of observations, Dec 18. 
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Figure 7. Four primary F814W-dropout candidates behind RCS2327. From left to right: Combined F435W + F814W optical image, 
F098M, F125W, F160W, combined F098M + F125W + F160W infrared image, IRAC [3.6X IRAC [4.5] , and rgb image where b = F098M, 
g = F125W, and r = F160W. Dropout IDs are shown in red on the rgb cutouts. All cutouts are 9" x 5.5'^ 


Table 4 

Dropout Candidates 


ID 

“J2000 

^J2000 

F160W 


•^phot 

D/i* 

logw{Agelyr) 

logio(M/MQ) 


RCS2327-1914 

351.878006 

-2.0875638 

25.6 

^•^^-0.24 

6.9 

r» qq + 0.03 
^•^^- 0.02 

Q r)<-)-|-0.02 

°-^^-0.03 

9.08lg;38 

2.2 

RCS2327-1282 

351.880689 

-2.0763742 

24.8 


7.1 

^•^^-0.05 

Q 7''7+0.02 

O- ' ' -0.04 

9 39+0 04 

1.0 

RCS2327-2068 

351.856179 

-2.0933244 

26.6 

10.42l?;Oi 

6.8 

^•^^-0.01 

Q 9 C + O.OI 
0.zo_q_24 

Q r)Q“l”0.32 
°-^°-0.25 

4.9 

RCS2327-1412 

351.856215 

-2.0785714 

27.0 

10.07+g;I2 

6.8 

u.uz^_o 002 

8.051^6 

Q ro + O.OO 
5.t)5_o 71 

7.7 


Note. — Properties of the four primary F814W-dropout candidates. The column F160W lists is the apparent magnitude in that 
filter. Photometric redshifts, z-^Yiot were obtained from the Le Phare Zml keyword. The full photometric redshift distributions 
are shown in Fig, fsl The characteristic luminosity, L*, was determined from the luminosity function at 2 = 6.8 determined by 
|Bouwens et al.|2Ul* log\{){M/MQy) represents the “demagnified” stellar mass. The final column represents the median rest frame 
XycTequIvaTentwTdfn noise (Icr). All errors reflect 68% confidence. 
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Figure 8. Spectral energy distributions (top panels) and photometric redshift distributions (bottom panels) for the four primary F814W- 
dropout candidates found behind RCS2327. Shown are fits to the photometric data points using HST data only (dotted blue line) and 
using HST+ IRAC [3.6] and [4.5] data (solid red line). 
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Figure 9. Line-like feature at AoUg = 10245 A in 3 hr combined 
MOSFIRE Y-band spectrum of dropout candidate RCS2327 - 843. 
If the feature is Lyo, which we find unlikely, the galaxy is at 
z = 7.42. The red horizontal lines indicate the expected verti¬ 
cal position in the slit from the broadband filter. The feature is 
well separated from nearby sky line residuals. The portion of the 
spectrum shown is 9" in the vertical direction. 


Furthermore, the photometric redshift distribution rules 
out z = 7.42 at greater than 98% confidence, even if we 
subtract off the line flux observed from the flux mea¬ 
sured in F098M in HST before computing the photomet¬ 
ric redshift. The photometric redshift distribution peaks 
at z = 6.7 and similarly prohibits [OIIJAA3726-3729 at 
z = 1.74, the case in which the photometric break is the 
Balmer break instead of the Lyman break. The night 
in which the feature appears had the worst seeing (Ta¬ 
ble [^. We therefore conclude that it is very unlikely 
that the feature is an emission line from RCS2327 - 843. 
RCS2327 - 843 is well isolated in HST data; the closest 
object that we find is at a distance greater than 3" from 
the position of the emission line, whereas the seeing was 
only ~ 1.20" during observations that night. The feature 
is therefore unlikely to be an emission line produced by 
a neighboring object. 

In addition to visually searching through all candi¬ 
date dropout spectra, we conduct an automatic search 
by extracting ID signal-to-noise spectra from the 2- 
dimensional reduced spectra. With a, S/N = 4 detection 
threshold in each separate dropout spectrum, 10 features 
are detected in the spectra of the four primary dropouts. 
None of these features pass tests 1-6, however. 

We calculate Lya flux limits in Y-band following 
the flux calibration and 1-dimensional extraction proce¬ 
dures described in Section Not accounting for mag¬ 
nification, we measure median 5a Lya flux limits of 
0.5 — 0.6 X 10~^^ erg s~^cm~^ , in ag reement with the 
limits observed by Treu et ah] (|2013 ), who took obser¬ 


vations with similar seeing and exposure times. Af¬ 
ter accounting for magnification, we achieve limits of 
5—14x 10“^® erg s“^cm“^. Rest frame EW limits are ob¬ 
tained from the flux limits relative to the HST continuum 
measured in F098M (Table|^ . Since EW is a relative flux 
measurement, it is unaffected by lensing magnification. 
We measure median rest frame EW limits of 1 — 8 (Icr) 
for the primary four dropout candidates. 

The fact that we do not detect significant Lya emission 
is unsurprising given our small sample size and the small 


et al.||2012 

Treu et al.||2012 

Bradac et al.||2012l h'inkel- 

stein et al. 

|2013||Treu et al.| 

20131 IPentericci et al.||20141 


Schenker et al.||2014). This is in contras t to a Lya frac- 
tion of ~ 5U7o at z ~ 6 ( Stark et al. 12010 1. The sharp de¬ 
cline in the Lya fraction from z = 6 to z = 7 is typically 
attributed to the attenuation of Lya due t o the rising 
abun danc e of neutral hydrogen in the IGM (|Treu et ah' 


2013 and Schenker et al. 2014). Mesinger et al.| (SOTS]) 
show with simulations of the IGM during the reioniza- 
tion process that such a sharp decline is unlikely due 
to IGM attenuation alone. Alterna tive explanations in- 
clude evolution in galaxy properties (Dijkstra e t al. 


or a rising abu ndance of absorption systems (Holton &: 


Haehnelt |2013 ). 


Stark et al. (2014) suggest using rest-frame UV lines, 


such as CHIJA1909A for redshift confirmation in place 
of Lya. While typically a weaker line than Lya, 
CHI]A1909A is not resonantly scattered by neutral hy¬ 
drogen. CHI] is of particular interest in our sample be¬ 
cause we probe the sub-L* population, from which 90% 
of the total UV luminosity density comes, at least at 
z — 2 dReddy fc Steidelj|26o 9[ |Oesch et al.||2010[ |Alavi 
et al.|2(H4[ [Stark et al.|2014[). \Ve conduct similar visual 


and automatic searches for the CHI] A1909A doublet in 
our H-band spectra of the dropout candidates. We run 
tests (1-6) on potential emission lines in H-band as we 
did in Y-band. Again, we see no convincing emission 
lines. 


5.1. Contamination 

M, L, and T dwarfs have similar colors to z ~ 6—7 star- 
burst galaxies (Fig. [^. While RCS2327 is well removed 
from the galactic plane {b ~ —58), we nonetheless inves¬ 
tigate M, L, and T dwarfs as a source of contamination 
in detail. We inspect each of the dropout candidates in 
HST to determine whether it is resolved. All of the can¬ 
didates are clearly resolved except RCS2327-1282, the 
brightest object in our sample. We attempt to fit its 
photometry to M, L and T dwarf spectra using the SpeX 
prism libraries^ We are unable to find a good fit to any 
of the dwarf SEDs. The SpeX libraries only extend out 
to the near-IR, but RCS2327-1282 is detected at > lOcr 
in both IRAC [3.6] and [4.5]. At these wavelength s, we 
rely on a s ample of M, L and T dwarfs studied by |Pat- 
ten et al. (2006) for the comparison. We compare the 
F'125W - [T^nd [3.6] - [4.5] colors of RCS2327-1282 
with the colors of the dwarfs, finding consistency with 
their T dwarfs only. Therefore, we cannot rule out the 
possibility that RCS2327-1282 is a T dwarf contaminant. 

Another source of contamination in our sample is the 
population of elliptical galaxies at z ~ 1. The addition 
of IRAC [3.6] and [4.5] to HST data can help lift the 
degeneracy between z < 2 galaxies and z > 6 starburst 
galaxies. Fig. illustrates this. In particular for can¬ 
didates RCS2327 - 1914 and RCS2327 - 1412, for which 
we obtained IRAC limits only, high redshift solutions are 
much more favored once the IRAC data are included. 

6 . CONCLUSION 

Our new gravitational lens model indicates that RCS2 
J232727.6-020437 is an excellent cosmic telescope, as ef¬ 
ficient as the Frontier Fields galaxy clusters. We discover 
16 new multiple images consisting of 6 new systems. We 

http://pono.ucsd.edu/-adain/browndwarfs/spexprism/ 































































13 


identify a strong emission line at 15436 A in the MOS- 
FIRE H-band spectrum of one arc, which is likely the 
[OIIIJA5007 emission line at z = Ao83. The spectro¬ 
scopic redshift suggested by the emission line is in concert 
with the tight constraints from the photometric redshift 
of the multiple image system. 

We find four highly magnified sub-L* F814W-dropout 
candidates behind the cluster. Measuring their photo¬ 
metric redshifts, we find that they are likely at z ~ 7. We 
also find four F814W-dropout candidates that fulfill less 
strict selection criteria. We search for Lya emission from 
all eight candidates using Keck/MOSFIRE. We identify 
a feature at 3.6a at 10245 A in the Y-band spectrum 
of one the secondary dropout candidates, but conclude 
that it is unlikely a true emission line. We see no con¬ 
vincing emission lines in the four primary dropout can¬ 
didates, reaching down to median 5cr Lya flux limits of 
5 — 14 X 10“^® erg s“^cm“^. We also measure median 
Icr Lya rest frame EW limits of 1 — 8. Because Lya 
emission is affected differently by the presence of neutral 
hydrogen in the IGM in different reionization scenarios, 
the Lya EW distribution at z > 7 can be used to distin¬ 
guish between models of reionization. The EW limits we 
measured from our MOSFIRE follow-up will be used as 
part of a larger sample obtained from dropout searches 
behind all SURFSUP clusters to help better understand 
the nature of reionization. 
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